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In this paper, we explore an idea of having Newton's constant change its value depending on the 
curvature scale involved. Such modification leads to a particular scalar-tensor gravity theory, with 
the Lagrangian derived from renormalization group (RG) flow arguments. Several of the well-known 
f[R) modified gravity models have remarkably simple description in terms of the infrared renor- 
malization group, but not the "designer" types in general. We find that de Sitter-like accelerated 
expansion can be generated even in the absence of cosmological constant term, entirely due to run- 
ning of the Newton's constant. In hopes of tackling the problem of cosmological constant's smallness, 
we explore the flows which are capable of generating exponential hierarchy between infrared and 
ultraviolet scales, and investigate cosmological evolution in the models thus derived. 
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I. INTRODUCTION 

It has long been hoped that quantum theory of gravity, 
at least in some limit, allows description in terms of an 
effective field theory [1] [2]. The usual Einstein- Hilbert 
action is merely the first two terms of an effective action 

/|E^'"'V"HA)7e(") + ...|y^d4x, (1) 

expanded in local n-th order curvature invariants 72. 
which explicitly depends on the ultraviolet cutoff scale A. 
Should the effective theory of gravity be asymptotically 
safe, this would allow a sensible ultraviolet-complete de- 
scription [3ii4j. However, recent observations of acceler- 
ated expansion of the Universe indicate the presence of a 
tiny but nonvanishing cosmological constant [5W8] , which 
presents technical hierarchy problem between (infrared) 
cosmological acceleration scale and (ultraviolet) Planck 
scale, which seems irreconcilable despite the many efforts 
put forward [9, 10 . 

Given enormous separation of cosmological and Planck 
scales, it would seem highly dubious that quantum grav- 
ity is somehow responsible for cosmological acceleration. 
Corrections and corresponding beta functions reviewed in 
[31 3] originate from ultraviolet degrees of freedom; they 
are Planck-suppressed and irrelevant for cosmology. But 
the language of effective field theory is universal, and 
should be applicable to infrared phenomena on cosmo- 
logical scales as well. To give cosmological acceleration 
a dynamical origin without invoking a new dark energy 
matter component, one would need a "modified" gravity, 
with new infrared degrees of freedom becoming active 
at cosmological scales, while remaining hidden in solar 
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system and local tests. A lot has been said on the sub- 
ject, with ,f{R) gravity models in particular receiving a 
lot of attention as of late [TT], yet the question of large 
hierarchy ever remains in these attempts. 

The main aim of this paper is to apply the idea of 
renormalization group, which has been very fruitful in 
high energy physics, to description of these new (and still 
unknown) infrared degrees of freedom. This obviously 
entails the shift of attention from ultraviolet cutoff in the 
effective action ([T]) to the lower limit of integration - the 
infrared cutoff , and on how things scale when it is varied. 
Of particular interest is the question whether it is even 
in principle possible to generate anomalously low scale of 
cosmological constant by specific running, perhaps in a 
way similar to dimensional transmutation phenomenon, 
which is what we will discuss here. 

Analyzing renormalization flows dependent on curva- 
ture scale, we find that de Sitter-like accelerated expan- 
sion can be generated even in the absence of cosmological 
constant term, entirely due to running of the Newton's 
constant. This presents a novel way to view the hierarchy 
problem, in contrast with previous studies jl2lll3) involv- 
ing running of A. Cosmological constant should be set to 
zero and protected by symmetry, while de Sitter asymp- 
tote is attributable to flow of Newton's constant, which 
we describe by renormalization group equations. Several 
historically important modifled gravity models of f{R) 
type |14H16) turn out to have a very simple description 
in this terms, but not the "designer" types, in general. 

We argue that a "soft" running, in which Newton's 
constant flows with quadratic beta function approach- 
ing a high-curvature fixed point, generates exponential 
hierarchy between infrared and ultraviolet scales, which 
could be exploited for cosmological model-building. The 
simplest model of this type, with beta function analo- 
gous to the one for QCD coupling constant, is somewhat 
problematic observationally (due to Newton's constant 
running to zero in high curvature limit), but is easy to 
analyze analytically. We present it here due to its sim- 
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plicity, and then discuss more realistic models with expo- 
nential hierarchy derived by this method, which did not 
appear in the literature before. 



Thus, Newton's constant in Einstein-Hilbert action be- 
comes a function of curvature, and the theory of gravity 
is promoted to /(i?)-type [11] with Lagrangian 



II. INFRARED RENORMALIZATION GROUP 
AND MODIFIED GRAVITY 

Let's explore the implications of a running Newton's 
constant with an energy scale set by curvature scalar R, 
which we will describe by introducing a dimensionless 
coupling constant a by 



a m 



pi 



(2) 



Note that presence of reduced Planck mass mp\ in the 
above definition means that gravitational coupling is 
compared to a fixed scale (namely the ultraviolet limit 
value) , which only makes sense if one is talking about its 
change with infrared cutoff, as we do here. We will as- 
sume that explicit cosmological constant term in effective 
action ([T]) strictly vanishes, and we will ignore higher or- 
der curvature corrections as we are primarily interested 
in infrared effects. If rcnormalization group flow is au- 
tonomous, the running of the dimensionless coupling with 
scale ^ = R/Rq is described by a beta-function flow 



(3) 



which is integrable for one-dimensional dynamical system 
da f d/i 



(4) 
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f(R)- 



(5) 



The crucial point here is that autonomous flow ([s]) is 
independent of Rq (which is merely a reference scale), 
and f{R) = R/a generated by it does not explicitly in- 
volve a tiny dimensionful scale Rm introduced by hand 
in most if not all f{R) models proposed so far [T5U20j . 
which replaces cosmological constant but of course still 
presents the same hierarchy problem. Instead, the hi- 
erarchy appears in particular solutions connecting local 
high-curvature environment (where we seem to live right 
on top of ultraviolet fixed point) and the cosmological ac- 
celeration attractor in the infrared limit. We will return 
to this point again in the next section, but for now let's 
explore f{R) models generated by such flows. 

For f{R) gravity to be perturbatively stable, one needs 
/' > to avoid ghosts {TP , and /" > to avoid tachyon 
instability [22,.23J. These conditions are readily rewritten 
in terms of flow equation quantities 



/' = ^, Rf" 



l-2^+/3'')4- (6) 



In addition, /(i?) gravity equations of motion ( 24|25 ) 
allow a vacuum (anti-) de Sitter solution if 



(7) 
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From ^ it follows that de Sitter solution is unstable if 
/3' < —3. If one wishes to have asymptotically safe theory 
with effective cosmological constant, the renormalization 
group flow should connect a + /3 — line in IR to a /3 = 
fixed point in UV, while avoiding forbidden region f3 > a. 
A number of such flows is shown in both panels of Fig. [T] 
both with stable and unstable de Sitter vacua. 

Interestingly enough, Einstein-Hilbert action with cos- 
mological constant, Einstein-Hilbert action with ultravi- 
olet corrections [M] , and original 1 /R proposal for in- 
frared corrections [151116] are all generated by very simple 
renormalization flows, as shown in Fig. [T^, all of which, 
jumping ahead of ourselves, only differ by a single co- 
efficient of proportionality (12). For example, requiring 
/' = 1 identically in we should recover general rela- 
tivity action. Indeed, integrating flow ^ with 



I3{a) — a~ o? 
leads to a rational dependence of a on curvature 



1' 



(8) 



(9) 



which puts back the constant term in the action 

i?-i?o- (10) 
All power-law corrections to Einstein-Hilbert action 



a 



f{R) = R 



1 + A 



R 
Rq 



are in fact generated by an autonomous flow 
/3(q;) = n a{a — 1). 



(11) 



(12) 



This could be checked by direct integration as we did 
above, but it is instructive to "reverse engineer" the flow 
( 12 1 from expression ( 11 ), to see how the scale i?o disap- 



pears from the flow equations. Writing a = R/ f{R) 

1 



1 + A^" ' 



(13) 



and formally taking logarithmic derivative /3 — fida/d^, 
we obtain the expression for /3 as a function of scale ^ 



nXfi" 



(1 + Am")2 ■ 



(14) 



This by itself is simply an identical rewrite of expres- 
sions for f{R) and its derivative /' in new variables, but 
if scale fi could be eliminated in favour of the coupling 
value a, the dynamical system would close and become 



autonomous. Inverting equation (13), we see that /3 in- 

(15) 



deed can be written as a function of a only 

/3 = -na^ia^^ - 1), 
which leads to the autonomous flow ([121 




FIG. 2: Not all f{R) models have a simple renormalization 
group flow interpretation - e.g. Starobinsky model jlSj shown. 

Even in one-dimensional dynamical systems of the type 
we consider, this inversion is not guaranteed to happen. 
Thus, not all f{R) models allow simple interpretation in 
terms of renormalization group flow. As a representative 
example, let's consider Starobinsky model [T5^ 



f{R) = R 



l + -((l + /i^)-"-l) 



R 
Rq 



(16) 



+ -((1 + m')-"-i), (17) 



This function corresponds to a non-invertible relationship 
between coupling strength a and curvature scale fj, 

-1 _ fi R) ^ A 
a — — — — 

R ^ 

which leads to a multi-valued function /3(a) as shown in 
Fig. [2] (for particular values of n = 1 and A = 2). Follow- 
ing the flow along cosmologically relevant branch (shown 
in red) from de Sitter solution toward increasing curva- 
ture, one jumps through R — > ±00 at a = 1, switching 
over to the negative curvature branch (shown in blue), 
and then passing through i? = at a = 1 again to finally 
reach small positive curvature branch (shown in green). 



III. IMPLICATIONS FOR HIERARCHY 
PROBLEM 

General behavior of modified gravity theories described 
by renormalization group fiow ([3| in high-curvature limit 
is determined by an asymptotic expansion of the flow 
near the ultraviolet fixed point a* 



/3(a) 



00 

E 

k=l 



Ckia - a*)'' 



(18) 
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FIG. 3: Running /3 as a function of scale ^ = R/Ro for linear 
(green) and quadratic (black) approach to UV fixed point. 



In the models we discussed so far, this series starts from 
a linear term (as shown in Fig. [T^), and integrating the 
flow leads to a power-law relationships between scale /i, 
coupling a, and running /3. With respect to cosmological 
curvature scale, there is little difference between local en- 
vironment on Earth and Planck scale; both are infinites- 
imally close to UV fixed point, and there are no obvious 
reasons why the scale separation should be so large. 

But if approach to UV fixed point a* is slower than 
linear (i.e. ci — 0) instead (as shown in Fig.[lj3), integra- 
tion of fiow equations leads to a logarithmic dependence 
of coupling a and running /3 on the curvature scale /i, 
which means that reasonable separations of the points on 
the flow diagram l3{a) can become exponentially large in 
physical curvature. Placing an upper limit on the run- 
ning /3 of Newton's constant in near-Earth environment 
(for example, by null deviation from Newton's force law 
in torsional pendulum experiments [24] ) . we then sud- 
denly have something to say about how far we should be 
removed from the infrared scale so as not to detect any 
running. 

This argument is illustrated by Fig. [3] which shows 
running /3 = da/dln/i as a function of curvature scale fi 
for two flows, /? = —a{a — 1) (in green) and /3 = —a^ (in 
black). For example, placing an (arbitrary) upper limit 
on < 10^'' does not tell us much about the first flow, 
yet requires a separation of at least 30 orders of mag- 
nitude in curvature between infrared and measurement 
scales for the second flow. In other words, deviations 
from Einstein gravity depend on how close we are to UV 
flxed point in the lab measurements, thus the smallness 
of cosmological constant is tied down to the observation 
that f{R) gravity reproduces Einstein theory! 

This simple thought could be used to generate small 
numbers involved in the hierarchy problem in cosmolog- 
ical model-building; the only thing we need to change is 
to make beta function quadratic in UV. To illustrate the 
general idea, let's consider the simplest model of such 
kind, namely the one where Newton's constant runs to 
zero with quadratic beta function in UV limit (which is 
similar to running of the QCD coupling constant) 



(19) 



Integrating the flow equation ( 19 ) leads to an effective 



action for gravity with 

fiR) 



R f R 

— 1-}- ao In — 
"0 V -f^o 



(20) 



Here ao is the coupling measured at curvature scale Rq^ 
and its value also determines how fast the Newton's con- 
stant runs at that scale. Note that i?o is not really a 
parameter in the action, but merely the scale at which 
the action is evaluated. Coupling a can be renormalized 
to another arbitrary scale i?i by redefining 



Ql = 



1 



«0 

ao In ^ 



(21) 



-Ro 



without any change to the Lagrangian itself. Corrections 
of the form ( 20 ) already appeared in [25l ESj , although 



full implications for hierarchy were not realized then. 



The obvious problem with the simplistic model ( 19 ) is 



that Newton's constant actually does run to zero, which 
of course is kind of different from the expected behaviour 
in large curvature limit. This drawback is easily rectified 
if one considers Newton's constant running to a constant 



/3(a) = ~n{a — a^Y 



(22) 



Integrating this flow leads to a rational function of loga- 
rithm of curvature of the form 



f{R) = 



6 In/i 



d In /i 



R 
Rq 



(23) 



As before, constants in the action can be redefined by 
changing the reference scale Rq, but the group action is 
slightly more complicated. Running of /3 with curvature 
scale is essentially the same as in simpler model (19 1, 
and is not shown in Fig. [3] to avoid crowding the plot. 
While more realistic, this model is more cumbersome to 
analyze analytically. As far as we know, it has not been 
discussed before. 



IV. COSMOLOGICAL EVOLUTION 

Now let's discuss cosmological evolution of the models 
we introduced in the previous section. Variation of the 
Lagrangian ([s]) with respect to the metric g^i^ yields field 
equations of motion in f{R) gravity 



^-,2r^„ (24) 



which involve higher derivatives of the metric, as usual 
in f{R) gravity. Taking a trace of the above equation 



□/' = J (2/ - f'R) 



Pi 3' 



(25) 



and introducing a scalar degree of freedom (f> = f — 2, 
equations ( 24|25 ) are cast into a scalar-tensor theory form 



G 



U4> = y (0) - J". 



(26) 




FIG. 4: Effective potential for scalar field 



Here Q^^ is the effective stress-energy tensor of geometric 
degree of freedom (jj (sometimes referred to as scalaron) 



Qm. = -(1 



G 



□0-3P(0) g^,, (27) 



with V'((j)) and P(4>) defined by equations (30|31 1, and 



J" 



T/3 is a matter force term for field 0. As 



Einstein tensor G^^ can be "traded" between left and 
right sides of equation, definition of Q^^ is not unique. 

The above applies for a general f{R) gravity theory. 
Equations for our model ( 20 ) become quite simple. Scalar 
degree of freedom cf) is just a logarithm of the curvature 



In 



R 

Rq 



1 - 
ao 



In 



R 



(28) 



where the constant in the definition is chosen so that 
(|>^. = solution corresponds to a de Sitter vacuum with 



i?* = 4A = i?o exp 



ao 



(29) 



Note that if the running parameter is small (ao <Si 1), 
effective cosmological constant A is exponentially sup- 
pressed compared to the reference scale i?o! In terms of 
A, the (/)-dependent terms in equations of motion ( 26|27 ) 
become 



P(0) 



(/ - f'R) = 
(2/ - fR) 



Ae"* 



(30) 
(31) 



Effective potential F(0) of the scalar degree of freedom 
is readily integrated and has a simple analytic form 



(32) 



The potential has a single de Sitter minimum at = 0, 
flat small curvature asymptotic (<^ — )■ —oo as i? = 0), an 
exponential potential wall for large curvature (0 — !■ +oo), 
as plotted in Figure |4] The matter force term drives the 
equilibrium field value V'{(j)cq) — up the potential wall 



= W 



3T 
4 A 



(33) 




FIG. 5: Poincare projection of dynamical system (351 



with explicit equilibrium position given in terms of Lam- 
bert M^-function. The field cj) is light near de Sitter min- 
imum, but for matter densities substantially higher than 
effective cosmological constant T/A ^ 1, field (p becomes 
heavy, and one expects the usual chameleon mechanism 
to "freeze" the scalar degree of freedom [27l [23 . 

Cosmological evolution of a flat Friedmann- Robertson- 
Walker universe with a metric 



-dt^ +a^{t) dx^ 



(34) 



in our model can be completely analyzed using dynamical 
systems techniques. The phase space is four-dimensional, 
with variables {0, tt, a, H} and dynamical equations 



3i7vr + y'(0) = mpt^, 



H=^- 2H\ (35) 
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supplemented with definitions (j) — n and d — aH. The 
evolution is subject to a constraint (which is analogous 
to Friedmann equation in the usual cosmology) 



+ 2) + i?7r + P(</)) ==m 



-2 P 
pi 3: 



(36) 



and hence is restricted to a three-dimensional surface in a 
phase space (which has a rather complicated shape). The 
clearest view of the trajectories is allowed by projection 
on a; = {(/), tt} plane followed by Poincare compactifica- 
tion y = x/{a'^ + x^)^ , as shown in Figure [s] Detailed 
analysis of the dynamics will be presented elsewhere, here 
we will just highlight the main features. The only finite 
critical point of dynamical system (351 is de Sitter at- 



tractor, which all cosmological solutions approach in the 



6 



future. Among all trajectories, the most important is the 
one tracking the diluting matter energy density 



'pi -? 



(37) 



shown in blue in Figure [5] In the past, it closely follows 
the minimum of the effective potential ( 33 1 due to matter 



W 



AI,0 



4A 



(38) 



until the scalar field (p becomes light and "releases" , ap- 
proaching de Sitter minimum kinematically. Tracing this 
trajectory to the present day value of constraint equation 
), shown as thick black line in Figure [5] gives a cos- 
ogical constraint on values of SIm, and A, similar 
to closure relationship in flat ACDM model. For famil- 
iarity's sake, one can introduce effective parameters 



(36 



no = 



PQ 



3m2ji72 ' 



WQ 



Pq 
PQ 



(39) 



describing scalar field contribution ( |27[ ) by extracting 

PQ^imlyH^-p, PQ ^ ml.iH^ ^ R/3) - p (40) 

from evolution history. As mentioned earlier, these are 
not unique, so for comparison with observations it is 
much better to directly calculate deviation of distance 
modulus from fiducial ACDM cosmology 



A/i — 5 log 



10 



"■L,f{R) 
di, ACDM 



1 



da 



(41) 



Assuming flM.o — 0.278 and matter-radiation equality at 
redshift z^q — 3250 [5], flatness prior fixes ^q q = 0.397 
at present for the tracking solution. Evolution with red- 
shift is shown in the left panel of Figure [6j While it 
is clear that expansion history in our model is substan- 
tially different from standard ACDM, deviation of dis- 
tance modulus up to redshift z = 2 is only about 1/20 
magnitude, which is not constrained by the present SN 
observations [7]. Also notable is the fact that effective 
equation of state of dark energy in this model can appear 
to cross phantom divide {w ~ ^1), and even develop a 
pole! 



V. DISCUSSION 



flow The bad news is that while the recent cosmolog- 
ical expansion history is plausibly reproduced, the New- 



ton's constant changes substantially (33 1 between cos- 



mological, galactic, and near-Earth environments in the 
simplest realization discussed here. Thus, although the 
scalar degree of freedom is not light inside matter, this 
model might have difhculty with gravity tests that probe 
absolute value of Newton's constant, for example con- 
straints on expansion during nucleosynthesis epoch. 

The difficulties with radically changing Newton's con- 
stant can be avoided by having a quadratic beta func- 
tion flow to a finite UV fixed point instead (22 1, which 
leads to models of the type (23). These models share 



similar features with respect to running and hierarchy as 
the simplest realization we discussed in detail above, but 
Newton's constant asymptotes to a finite value in the re- 
gions of high curvature, as in Einstein's gravity. Taking, 
for dcfiniteness, a* = 1/4 and k — 4/3, leading to 



1 



1< 



1 



i?, 



(42) 



one can repeat the analysis of the previous section. While 
algebraically more complicated, the story goes pretty 
much along the same lines. Effective potential V{(\>) 
still has a single minimum and an infinite potential wall 
at high curvature, but develops a turning point at the 
low curvature limit. For sensible cosmology, this turning 
point must be below asymptotic de Sitter curvature in 
value, which places certain condition on the fiow (22 1. 



The constants in (42) are chosen so that this is the case. 



otherwise no effort was made to tune them to "special" 
values or to fit ACDM fiducial model. With these values, 
the tracker solution is parametrized by JIq = 0.722, and 
is shown in the right panel of Fig. |6] 

The deviation of distance modulus from fiducial 
ACDM model is about 1/10 of magnitude out to red- 
shift 2 = 2, and f2 parameters plot looks much more con- 
ventional now. The tell-tale sign of this model is much 
slower dilution of fig at high redshift than usual, with 
il^f never quite reaching 1 during matter domination in 
this example. This will surely affect the structure for- 
mation, and should give a way to test or rule out the 
model. While parameter space is larger and still remains 
to be completely explored, this "soft" approach to Ein- 
stein's gravity is due to logarithmic dependence and is 
characteristic for all the models we presented here. 

Severity of bounds placed by solar system tests [T7| [25] 
and large scale structure growth [5DH5^ will also need to 
be investigated. 



The main result of this paper is that running Newton's 
constant can cause accelerated expansion of the universe. 
The good news is that the simplest model of this type 
( 19|20 ) is predictive and has essentially the same number 
of parameters as standard ACDM, with exponential hi- 
erarchy between effective cosmological constant ( 29 ) and 



UV scales naturally generated by renormalization group 



Acknowledgments 

This work was supported by the Natural Sciences and Engi- 
neering Research Council of Canada under Discovery Grants 
program. JQG was supported by the Pacific Century Gradu- 
ate Scholarship of the Province of British Columbia and the 



7 




30 
20 
10 

n 

-10 
-20 
-30 

0.5 



-0.5 













10-' 


10" 


10' 


10" 


10= 




























WQ - 


1/3 






10-1 


10» 


.-■'I'o' 10^ lO' 


10" 


10= 






■ \ 1 +z 










^ approximation V 


(0oq) 















(a) (l + ln^)i? 



1.6 
1.0 

^ttac 0.6 





present 




10-' 


1 


y 10' 


10^ 


10' 


10" 


10= 



1 +z 




0.5 



-0.5 



10-1 


10" 


lO' 10^ lO' 

1 +z 


10" 

WQ ~ 


105 

1/3 




10-' 


10" / 


lal-— — TT)^^ 10' 
<^ 1 +z 


10" 


10= 


WQ 




^ approximation V' 


(<^oq) 





(b) /(i?) 



1 + In■ 



l + 1 In ^ 
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